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Atmospheric Muon and Neutrino Flux From 3-Dimensional Simulation 
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The atmospheric muon and neutrino flux have been simulated using the same approach which 
successfully accounted for the recent secondary proton, electron and positron flux measurements in 
orbit by the AMS experiment. For the muon flux, a good agreement is obtained with the CAPRICE 
and HEAT data for altitudes ranging from sea level up to about 38 km. The general features of 
the calculated atmospheric neutrino flux are reported and discussed. The flux obtained at the 
Super-Kamiokande experiment location are reported and compared with other calculations. Eor 
low neutrino energies the flux obtained is significantly smaller than that used in the data analysis of 
underground experiment. The simulation results for the SOUDAN experiment site are also reported. 


I. INTRODUCTION 

With the rapidly increasing amount and statistical significance of the data collected by underground neutrino 
detectors Hi, the precise calculation of the atmospheric neutrino flux is highly desirable, since the interpretation 
of the data relies on the calculated flux as a fundamental input. In the contained events analysis, the ratio of 
ratios, namely the ratio of the observed number of muon events to the number of observed electron events with the 
corresponding ratio of the same numbers from simulation calculations, is not sensitive to the absolute neutrino flux. 
However, in the analysis of upward going muon events, the absolute flux need to be known accurately in addition 
to the detailed features of the flux such as the zenithal and azimuthal angle distribution Hi The availability of 
calculations of proven validity would discard all possible doubt on the conclusions. 

In the past 20 years, various methods including analytical, semianalytical [3, kinematical H and Monte Carlo 
techniques Hi, have been applied to calculate the atmospheric neutrino flux iim Ilin i2j. However, the 
accuracy of the results was limited by the accuracy on both the cosmic ray (CR) abundance, and the production cross 
section of the parent particle (tt^, in the collisions of CRs with atmospheric nuclei over the relevant 

energy range these calculations. These two quantities being basic inputs to the calculations. The 

analysis and the comparisons made in indicated that the particular choice of the primary spectra and the strong 
interaction model, as well as the earth geomagnetic field could result in significant differences in the calculated neutrino 
flux. 

The new AMS M and BESS m measurements of the CR proton flux consistently agree to within 5 %. The 
values reported are lower than those used in the previous works by about 20 to 25 % 0|. Besides, some analysis 
showed that the pion production cross sections given by the general purpose softwares based on theoretical models, 
may largely depart from experimental data. 

Although the values of the flavor ratio reported in previous works are very close to each other, the absolute values 
of the flux obtained in these works fail to be consistent with each other within satisfying precision In particular, 
these results seem not to predict quantitatively the observed East-West (EW) asymmetry of the neutrino induced 
events which is expected to originate from the EW asymmetry of the geomagnetic cut off (GC) on the incoming 

CR - mainly protons and helium particles here - momentum, due to the earth magnetic field. Any calculation, to 
be reliable, must reproduce the experimental EW asymmetry since flavor oscillations do not change the direction of 
motion. This requirement can be taken as a testing ground for the various approaches of the neutrino flux calculations. 

With the occurrence of the new measurements of the primary CR flux, together with the continuous improvement 
of the cross section calculations, most approaches have been updated [T^ , and new calculations have 

been proposed In parallel, experimental groups also made greater efforts aiming to understand 

the neutrino production process in the atmosphere with the measurements on the muon flux at various 

altitudes, together with those of the primary proton and helium flux. These new data provided a sounder testing 
ground for the reliability of the numerical approaches to the particle production and to the particle dynamics and 
kinematics in the atmosphere and in the earth geomagnetic field. 
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Particle type 

P 

He 

CNO 

Ne-S 

Fe 

< 4>{<z>) > 

106 

16 

1.1 

0.34 

0.08 

A 

V 

0.3 

0.63 

1.23 

1.53 

1.86 

fractional tt flux 

0.72 

0.22 

0.032 

0.011 

0.0035 


TABLE I: Mean flux 0 ( 2 :) {w? s sr) tt multiplicity </i(7r) >; of the components of the CR flux, and their product normalized 
to 1, for the conditions described in the text. 


The work reported here is an extension of a research program which primary motivations were to interpret the 
AMSOl measurements of the charged particle flux in near earth orbit Ullll. To this purpose, an event generator 
describing the CR induced cascade in the atmosphere, particle propagation in geomagnetic field, and interactions 
with the medium, has been developed and successfully used to reproduce the proton |^, electron-positron and 
helium 3 1^ flux data measured by AMS and their relevant dependence on the geomagnetic coordinates. Since the 
generator of the program was basically the same as needed to generate the muon and neutrino flux, the code could 
be rather straightforwardly extended to describe the latter and to address the important issue of the atmospheric 
neutrino flux. It should be emphasized that the present approach allows an evaluation of the neutrino flux induced by 
both the primary CR flux and the atmospheric secondary flux, which have been tested sensitively and independently 
on other observables, conferring to the results a character of sound reliability. The above issues will be discussed in 
details below. This approach is also being successfully applied to the calculation of the atmospheric secondary flux of 
antinucleons 

The article reports on the calculation of the atmospheric muon and neutrino flux. The paper is organized as follows. 
In section m the method and models used in the calculations are introduced. Section irm is devoted to the results on 
the muon flux. In section HVl the properties of the calculated neutrino flux averaged on the whole terrestrial sphere are 
discussed, while the properties of the local neutrino flux calculated at the location of the Super-Kamiokande detector 
are discussed in section 0 In the following section ED the neutrino flux at the Soudan detector location is discussed 
along the same lines. Summary and conclusions are given in section ivni 


II. METHOD AND MODELS 

The method of calculation and the models used in the simulation program are described in mil- The calculation 
proceeds by means of a full 3D-simulation program. The main features of this approach are listed below for the 
reader’s convenience, with emphasis on some particular points relevant to the neutrino and muon flux. 

1. CR generation: Incident Cosmic Rays are generated on a virtual sphere chosen at some distance from the earth 
as discussed below. Events are generated randomly and uniformly on this sphere. In order to get an isotropic 
flux at any point inside the volume of the virtual sphere, the differential element of the zenith angle distribution 
of the particle directi on g enerated on the sphere must be proportional to cos Oz d{cosOz)^ Oz being the zenithal 
angle of the particle jllL (see also Q). The particle {A and Z) and its momentum are then generated 
according to the CR abundances and spectra discussed below. Generating random events on a virtual sphere 
far enough from the earth, i.e., at a distance where the geomagnetic field effects are negligible, so that the 
CR flux is isotropic, would imply a very large distance from earth and thus the generation of a tremendous 
number of events, most of them being useless for the simulation purpose since they would not reach the earth- 
atmosphere system. Instead of this time consuming direct method, the generation sphere is chosen close to 
earth, at a 2000 km altitude. The geomagnetic cut-off is then applied by back-tracing the particle trajectory in 
the geomagnetic field. In this procedure, only those particles reaching a backtracing distance of 10 Earth radii 
are kept in the sample. Flux conservation along any allowed particle path in the geomagnetic field is ensured 
by Liouville’s theorem application 

2. CR abundances: For the incident CR proton and helium flux, functional forms fitted to the 1998 AMS mea¬ 
surements were used EHi . For other periods of the solar cycle than those of the measurements, the incident 
cosmic flux are corrected for the different solar modulation effects using a simple force law approximation |4n| . 
The main components of the incident CR flux are usually divided in the following most abundant groups of ele¬ 
ments: p. He, CNO, Ne — S, and Fe. The contribution of each group to the neutrino flux scales roughly with the 
product (/){< z >)jii{7r) with (l){<z >) the CR abundance of the element or group of elements with mean electric 
charge < 2 : >, and /i(7r) the average multiplicities of pious produced in the CR collisions with atmospheric nuclei. 
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This is approximate since secondary protons induced by the He component also contributing to this production 
are not taken into account in this product. Tabled shows the relative abundances of the estimated pion flux 
induced by the various components of the CR flux. They have been evaluated using the AMS measurements for 
p and He and those from for the heavier elements with the spectral index from The flux have been 
calculated for particle rigidities above a geomagnetic cutoff of 10 GV. The pion multiplicities were taken from 
^3,0 at 4.2 GeV/c per nucleon. The sum of the contributions from the A > 4 flux components amounts to 
about 6%, while they are 72% for p and 22% for He. This choice for the GC maximizes the He (and heavier 
elements) over p fraction, which would be smaller for a lower cutoff (the He fraction would be around 10% for 
no cutoff). The A > 4 components of the CR flux were not taken into account in the present calculations. They 
will be included in the further developments of the code. 

3. Particle propagation: Each particle is propagated in the geomagnetic field and interacts with nuclei of the local 
atmospheric density. The specific ionization energy loss is computed for each step along the trajectory. The 
model used to describe the atmospheric density was taken from references 

Every secondary particles are processed the same way as their parent particle, leading to the generation of 
atmospheric cascades. 

4. Secondary particle production: Nucleons, pious and kaons are produced with their respective cross sections. 

Eor these particles the Kalinovsky-Mokhov-Nikitin (KMN) parametrization of the inclusive cross sections 
was used. Eor pions the KMN parameters have been fit on a wide range of p 4- A tt^’^ data between 1.38 
and 400 GeV/c incident momenta Ullil- Leptons - electrons, muons, and neutrinos - are produced in the 
decay chains of mesons, mainly pions and kaons, hadronically produced on atmospheric nuclei A in processes 
of type: CR + + M ^ + + + See below for details. It must be emphasized here, 

that the constraint to reproduce simultaneously the lepton population measurements at various altitudes and 
geographical coordinates, is extremely strong, and that a successful result would ensure a high level of reliability 
of the theoretical grounds and of the numerical approach to the problem. 

5. p decay: Eor the decay of muons, the spectra of the products (z^, P, e^) are generated according to the Eermi 
theory. 

6. K decay: Eor kaon decay, the Dalitz plot distribution given in ^3 used. 

7. p polarization effects: Muon polarization was taken into account for pion decay and for the 
channel Umiii. 

8. CR Energy range: The kinetic energy range of incident CRs covered in the simulation is [0.2, 2000] GeV, from 
around the pion production threshold up to a value where the corresponding neutrino flux produced becomes 
negligible. 

9. Processing: Each particle history, trajectory parameters and kinematics, is traced and recorded by the program. 
The event file was then analyzed separately to generate the various distributions of interest. 

10. z/ detection: Eor the neutrino flux calculation, a virtual detection sphere is defined at the Super-Kamiokande 
detector altitude. 

11. Elux Normalisation The differential particle flux produced in the simulation (m“^ • s“^ • sr“^ • GeV~^) were 
calculated by dividing the number of particles detected by the (local) surface of the particle collection area, 
solid angle, energy bin size, and equivalent sampling time of the CR flux. The latter was obtained from the 
total event number generated in the simulation run(s) divided by the surface of the generation sphere times the 
integrated (cos Oz weighted) solid angle (tt) times the energy integrated flux. 

To increase the statistics for the energy regions of low particle flux, the primary spectra were separated into 
a set of energy intervals such as [0.2,50], [50,100], ..., [450,500], [500,600], ..., [900,1000], [1000,1250], ..., 
[1750, 2000] GeV, the simulation program was run separately for each primary energy bin, and appropriately 
renormalized later when the full sample was constituted. 

The corresponding exposure times ranged from 2.1 x 10“^^s for the lowest energy bin to 4.7 x 10“^^s for the 
highest energy bin, corresponding to a total of the order of 10^ events generated in the sample. 

These calculations as those reported previously mill 113 include no adjustable parameter. 
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U\JON FLUX AT VARIOUS ALTITUDES 
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FIG. 1: Left: Simulation results (histograms) for the negative muon flux at various altitudes in the atmosphere, compared to 
measurements (full circles and triangles), from sea level up to about 38 km. See text for details. Right: Ratio of simulated flux 
to the measurement for each data point. 


III. ATMOSPHERIC MUON SPECTRA 

As it was mentioned above, atmospheric muons are produced in the decay chain of pions and kaons produced in the 
collisions, i.e., in the same reaction chain as neutrinos. A successful reproduction of the muon flux in the atmosphere is 
thus a further support to the reliability of the neutrino flux calculated in the same framework, physical conditions and 
computational environment. This is true as well for the electron-positron flux. For the latter, the present approach 
has allowed to successfully reproduce the AMS measurements of the flux and of their latitude dependence, in near 
earth orbit recently (see for the in the atmosphere). The results on the muon flux are presented in this 
section. 

It must be emphasized first that the muon experimental observables, namely the energy dependence of the flux at 
various altitudes, are not available for neutrinos. Therefore they offer a test ground for flux calculations which is more 
sensitive to the various elements of the calculations than the neutrino observables discussed below. For this reason 
they deserve to be investigated with care. 


A. Negative muons 


Figure HI shows the calculated muon flux compared to the data measured by the CAPRICE and HEAT experiments 
at various altitudes (left), and the ratios of the simulated to measured values for each data point (right). The 
agreement between simulation results and data is quite good through the energy range investigated (0.5-50) GeV for 
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FIG. 2: CR proton kinetic energy distributions of the cumulative probability for producing downgoing muons in the energy 
bins of 0.2-0.4, 0.8-1.2, 2.0-4.0 and 8.0-12.0 GeV (from the left to the right), with zenith {Ozenith < 15°), for the altitudes 
indicated. 


(GeV) 

38 

26 

13 

Altitude (km) 

10 6 

1.27 

0.66 

0.36 

0.2-0.4 

8.3 

9.0 

14.0 

18.6 

35.6 

66.1 

63.6 

68.1 

0.8-1.2 

17.5 

16.6 

25.0 

30.8 

46.8 

75.9 

84.9 

82.0 

2.0-4.0 

33.0 

34.1 

46.9 

53.7 

74.7 

100.3 

102.9 

101.9 

8.0-12. 

93.2 

97.9 

127.2 

136.3 

142.3 

186.8 

189.5 

190.3 


TABLE II: Mean incident GR proton energy (GeV) producing near vertical {Ozenith < 15°) downgoing muons in the specified 
energy bins and altitudes. 


all altitudes, from about the sea level to 38 km. It is especially good over the region from 10km to 26km altitude 
where a large fraction of the neutrinos detected by underground detectors are produced (see below). A significant 
difference between simulation results and data is observed however at 38 km altitude for muon energies below 1 GeV. 
This might result from an underestimated production of low energy pions for CR incident energy below about 10 GeV, 
since the pion production cross section in this energy region is poorly constrained by the very little data available. 
Note however that relatively large uncertainties and variations of the measured values, are observed for the low energy 
muon data M- 

The energy distribution of the cumulated probability for the incident CR protons to produce downward-going muons 
close to zenith in different energy bins and for different altitudes, is illustrated in figure O The corresponding mean 
values of the primary proton energy are tabulated in table nn 

From figure El it is easy to see that the cumulative probability curves shift gradually toward the high energy region 
with the decreasing altitude, while the distances between the curves are compressed. This can be clearly seen as well 
from table im At 38 km altitude, the mean energies of primary protons producing 0.2-0.4 GeV muons and producing 
8.0-12.0 GeV muons differ by a factor of about 11, this factor becomes about 7 at 10km, and 3 at around sea level. 

Figure El table DU and the following figure El and figure IT^ show that, 

1. at the highest float altitude, the simulated muon flux is most sensitive to the pion production cross section used 
in the simulation and hence, the higher altitude muon flux data can thus be taken as a good test of this cross 
section. 
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FIG. 3: Flux ratios for muons induced by the primary CR (proton + helium) flux (solid line), and by atmospheric secondaries 
(dotted lines), respectively, to the total flux, for three detection altitudes. The total flux are those from figure Q] 


2. in contrast, at around sea level, the low energy muon flux should be more sensitive to the proton/neutron 
production cross section because on the average: 

(a) The low energy muon path length is roughly on the km scale, and the muons measured at sea level are 
thus produced at relatively low altitude. 

(b) Above 10 km, about 5 collisions have occurred. In these collisions, more secondary protons/neutrons are 
produced, that can contribute largely to muons parent meson production. 


B. Secondary proton contribution 

The ratios of the muon flux originating from primary CRs and from secondaries respectively, to the total flux, are 
shown on figure |31 The flux are nearly vertical and the corresponding normalization area are the same as those shown 
in figure Q for the three altitudes. 

This figure shows as it could be expected from simple collision rank considerations, that in the whole muon energy 
range, at the highest float altitude, more than 98% muon flux originate from primary CRs. At intermediate altitudes, 
the flux induced by secondaries decreases continuously with the energy from about 60% at 0.2 GeV down to about 
20% at 10 GeV, while at around sea level, the same trend is observed with a slower decrease with energy. 


C. Positive muons 

Figure 0 shows the calculated /i+ flux compared to the measured data from the CAPRICE and BESS experiments. 
The agreement with the BESS 99 data at mountain altitude is very good, but in the low (below 1.0 GeV) and high 
(beyond 10.0 GeV) energy range, the departure is obvious. This may result from the overestimation of the secondary 
proton/neutron production cross section at low energy and pion/kaon production cross section at high energy. 
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FIG. 4: Left: Simulation results (histograms) for the positive muon flux at various altitudes in the atmosphere, compared with 
measurements (full circles and triangles), from sea level up to 38 km, from refs|^ . 1 ^ . See text for details. Right: Ratio of 
simulated flux to the measurement for each data point. 


IV. PROPERTIES OF NEUTRINO SPECTRA 

In this section, the characteristics of the simulated neutrino spectra, averaged over the whole detection surface, for 
three bins of geomagnetic latitude and integrated over 4 tt solid angle, are described in details. The detection sphere 
is defined at 0.372 km, the altitude of the Super-Kamiokande detector. Some features having important implications 
for underground experiments, are emphasized. 


A. Absolute flux and flavor ratios 

The calculated flux are shown in figure [3 for three bins of geomagnetic latitude. The corresponding ratios of flux 
are shown on figure El where it is seen that the flux increases with the geomagnetic latitude in the low energy (< few 
GeV) range. At 0.1 GeV, the difference between the flux at low and high latitudes can be more than a factor of 2. 
However, this difference decreases with the increasing particle energy, vanishing beyond a few GeV. This reflects the 
effect of the geomagnetic field. Indeed, the mean primary particle energy associated to 1 GeV neutrino production, 
is about 60 GeV (see below), a value at which the effects of the geomagnetic cutoff (GC) are small. 

Figure El shows that the flavor ratios for low, middle and high latitudes are very similar to each other. The is 

flat through energy range, with its value near about 1.2. Vej^e decreases steadily from nearly 2 to about 1.2 - overlap 
with at high energy end, while {y^ + yn)/(yp. + z^e) rises up steeply above about 3 GeV. Large difference for 

ye/yp from the previous works is seen here (See |3illl| for comparisons). 
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FIG. 5: Neutrino flux in three bins of geomagnetic latitude 0-0.5, 0.5-1, and l-7r/2 radians, averaged over 4 tt solid angle. The 
error bars given in the figure and all the following figures are statistical unless specified otherwise. 
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FIG. 6: Individual vjv and sum (y^ + + ^e) flux ratios, as in figure|^ 


These features can be easily understood in terms of the productions of 7r+(iF+) and and of the muon 

kinematics. The ratio Vet^e is principally controlled by the ratio of 7r+(iF+) to 7r~{K~). The ratio {u^ + yii)j{ye + 
Pe) is sensitive to the muon decay kinematics. The variation of with energy is thus determined by both 

7r+/7r“(iF+/iF“) production ratios and the muon decay kinematics. Because of the Lorentz time dilation, more 
higher energy muons do not decay while crossing the atmosphere. Hence, the z^e(z^e) production (resulting from muon 
decay) is lesser than for z/^(z/^), since a has been produced in association with the y in the tt decay, resulting in a 
(z^/x + + ^e) ratio rising up with the neutrino energy. 

Here, two points should be noted. First, the z^e/^e ratio at low energy should be expected to be larger since the 
largest production ratio of tt^/tt” is at low incident CR energy (see for example). Second, the z/g/^e ratio should 
decrease to overlap with at high energy with their values approaching to 1.0. This is because at high incident 

primary energy, the production ratio of tt^/tt” and K~^fK~ approaches to 1.0 lEi, and most of the high energy 
muons from which the high energy z/g and Pg come will not decay before reaching the detection sphere. So, the z/^ 
flux is dominated by the 7r+ decay while the flux is dominant-ed by 7r“ decay. And hence, at very high energy, 

Ve/Pe ~ V^IPf, « 7r+(K+)/7r"(K“). 

Such characteristics can be seen clearly in figure [SI 
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FIG. 7: Left up and middle: Energy spectra of the CR Helium flux contribution (dashed lines) and total Kaon decay contribution 
(dotted lines), compared to the total flux (solid lines), for the v and z/, neutrino flavors from top to bottom as indicated in the 
figure respectively. Left bottom: Ratios of CR Helium induced + i7^) flux to total + i7^) flux (solid line), CR Helium 
induced {ue + ^e) flux to total (ye + Fe) flux (dashed line), kaon decay induced (v^ + i7^) flux to total (v^ + i7^) flux (dashdot 
line), and kaon decay induced (h>e + Fe) flux to total (ue + Fe) flux (dotted line). Right: Flavor ratios for the same three 
contributions (Total, CR He induced, and Kaon decay, from top to bottom respectively), averaged over 4 tt solid angle and 
integrated over the whole detection sphere. See text. 


B. Contribution from the incident CR He flux 

The calculated total neutrino flux is compared on left figure [7| with the contribution induced by the CR Helium 
(He) flux and with the total kaon decay contribution. The corresponding flavor ratios are also shown on the same 
figure (right panel), where it is seen that the contribution of the CR He flux to the total neutrino flux is about 10% 
through the whole energy range. 

This contribution is constrained by both the CR He spectral abundance and the He induced pion multiplicity 
in nuclear collisions. Around 20 GeV per nucleon incident energy, where the neutrino production around 1 GeV is 
maximum (see discussion below and [13): Ihe CR He flux is about 5% of the p flux, this fraction decreasing slowly 
with the increasing primary energy. Experimentally, the pion multiplicity increases with the mass of the projectile 
according to a known law [^ . For the a ^ C reaction at 4 GeV/nucleon, it is about 3 times the p C value at the 
same incident energy M- These numbers are in qualitative agreement with the 10% contribution to the neutrino 
flux found to arise from the He flux, assuming that the proton induced helium contribution is small. This result is 
also consistent with the previous work on the AMS lepton data by the authors M- 

The difference between the contributions of CR He to the (ue + Ue) and (u^ + flux ratios is very small as seen 
on lower left figure [3 For the flavor ratio of the He flux contributions shown on middle right figure [3 they have a 
similar energy dependence as the total flux, as it could be expected (top). 
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C. Contribution from secondary Kaons 

At low energies the neutrino flux mostly results from tt decay in account of the dominance of the pion production. 
The Kaon decay contribution increases however with the increasing incident energy and the K production cross 
section, and finally dominates the production for very high energy neutrinos. This is shown on left figure [7| where the 
kaon decay contribution to the total neutrino flux (bottom left) is seen to increase continuously from about 2-3% at 
0.1 GeV up to ~15% for (z/g + and ~20% for {y^ + above 20 GeV. This results is in agreement with 

The energy dependence of the flavor ratio for the kaon decay contribution to the neutrino flux is a little different 
from the averaged value (bottom right panel on figure[3). The ratio z^e/^e decreases slowly with the increasing energy 
from 1.7 at 0.1 GeV down to about 1 at 10 GeV, while increases continuously from 0.9 at 0.1 GeV up to about 

2.1 at 10 GeV. In contrast with the individual flavor ratios, the + yii)j{ye + ^e) ratio grows rapidly from 1.8 at 
0.1 GeV up to about 5 at 10 GeV. These differences principally result from the to K~ production ratio being 
larger than for 7r+ to 7r“ 

Another noticeable feature of the kaon decay contribution is that, below 0.4 GeV neutrino energy the ratio + 
+ ^e) is less than 2. This can be understood from the decay mode having always larger a branching 

ratio than the mode. In particular, for the K^^y channel is the dominant decay mode |4^ . 


D. Contribution of atmospheric secondaries to the neutrino flux 

The contribution of the atmospheric secondary flux is compared on figure |HI with that of the primary component. 
It is seen that, at 0.1 GeV, about 50 % of the neutrinos come from secondary proton/neutron induced reactions, this 
fraction decreasing to about 20% at 50 GeV. 

The ratio of the secondary contribution, z^e/^e is larger than for the primary contribution. This is because of the 
larger flux of low energy secondaries due to the larger tt^/tt” production ratio at low incident energy. 

The distribution of the rank in the cascade of the neutrino producing collision is plotted on figure EM The 
rank is defined as the number of collisions taking place before the neutrino is produced. It is 1 when the neutrino is 
produced from the incident cosmic ray interaction. The figure shows that, on the average, about 2.4 collisions occur 
for the parent pion or kaon to be produced. 


E. Zenith angular distribution and geomagnetic latitude dependence of the neutrino flux 

The zenith angle distributions are shown on figures 11171 ITTl and m\ for the geomagnetic latitude bins (rad) 0-0.5, 
0.5-1 and 1-|^, respectively, and for several energy bins between 0.1 and 10 GeV. 

The figures show that the zenith angle distributions for downward going neutrinos (cos ^ > 0) have about the same 
shape for all the energy bins, while for upward going neutrinos, they display quite different shapes for the low energy 
bins, showing: a) A clear maximum for low latitudes, b) A fiat shape at middle latitude, and c) A deep minimum at 
high latitudes. This can be understood by considering the following two particular cases: 1) For a detection point at 
the north pole, the (7r/4) zenith angle cos^ ^ — l/v^ points to particles coming from the equatorial region where the 
geomagnetic cutoff for CR particles is highest. The neutrino flux is then expected lower for this angle. This effect is 
more sensitive for the low energy region of the neutrino spectrum. 2) Conversely, for a detection point at the equator, 
the zenith angle cos^ ^ — \/2/2 points to particles coming from the poles (although for a narrow azimuthal angle 
region only) where the GC is minimum and hence the observed neutrino flux is expected larger around this angle. 

The zenith angle dependence of the flavor ratios is modulated clearly by the tt^/tt” production ratio, GC and the 
muon decay kinematics. Interestingly, in low energy bins, the z^e/^e dependence on the zenith angle displays a similar 
behavior as observed for the flux. In addition to its dependence on the GC, this feature is also related to the pion 
multiplicity. Because the production ratio for 7r+ to 7t~ is larger below 10 GeV , the incident CR flux within this 
energy range which contributes a large fraction of low energy neutrinos, is very sensitive to the GC, and hence the 
z^e/^e changes accordingly with the geomagnetic latitude and zenith angle. 

The VyilvyL distribution is fiat in low energy bins but behaves similarly as + ^e) at high energies. 

For both the downward and upward flux, it is expected that the ratio approaches asymptotically the I'k~ 

production cross section ratio at the limit of infinite muon lifetime. 

The (z/^ + ^^)/(z^e + ^e) ratio displays a clear dependence on the zenith angle, with a symmetric up-down distribution, 
increasing from a minimum around the horizontal direction to a maximum for vertical incidences, the larger the energy 
bin, the more dramatic being the variation. This is because muons coming in horizontal directions have long path 
lengths in the atmosphere and then enough time to decay, producing a z/^ — z/^ pair and a Uei^e) for each pion decay. 
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FIG. 8: Total neutrino flux, primary CR proton and helium directly induced neutrino flux, and the secondary induced neutrino 
flux, averaged over whole detection sphere and over Att solid angle. Top left: CR Primary (dashed lines, up: down: z/e), 

secondary (dotted lines, up: down: v^) contributions for and and total Ve and (solid lines, up: down: v^. Top 

right: CR Primary (dashed lines, up: z7^, down: z7e), secondary (dotted lines, up: z7^, down: Ve) contributions for z7e and z^, 
and total and (solid lines, up: down: v^. Middle: Ratios of Primary contributed flux to total flux and Secondary 

contributed flux to total flux, for v (left) and v (right). Dotted line: ratio of primary directly induced {i7^) to total 
(z7^). Solid line: ratio of primary directly induced Ue (j^e) to total Ue (Fe). Dashdot line: ratio of secondary induced (z7^) 

to total (C^)- Dashed line: ratio of secondary induced Ue (j^e) to total Ue (Fe). Bottom: Flux ratios z^e/z^e, and 

+ ^/^)/(z^e + ^e) for the Primary directly induced (left panel) and Secondary induced (right panel) components. 


The situation is different for muons moving downward or upward, since the flight path available is much shorter. In 
this case, a fraction of high energy muons do not decay before reaching the earth and then do not contribute to the 
z/e(^e) flux, wherefrom a larger (z/^ + + ^e) ratio. The larger the muon energy, the larger the effect. 


F. Azimuth angle distributions and geomagnetic latitude dependence 

The azimuth angle distributions of the neutrino flux and of the flavor ratios are shown on figures IT^ to ca in bins 
of neutrino energy (0.1-0.3, 0.3-1, 1-3.1, 3.1-10 GeV), and for three bins of geomagnetic latitude (in absolute value): 
0-0.5, 0.5-1, and l-7r/2 rad, respectively. 
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Number of colHsIons 


FIG. 9: Rank distribution of the neutrino producing collisions in the atmosphere. See text for details. 


It can be observed on these figures that the azimuth angle distributions display a significant dependence on the 
latitude considered. For low geomagnetic latitudes (figure E|), the distributions are dominated by two peaks located 
symmetrically at the East (E) and West (W) azimuth angles, but with different heights. These latter features are 
clearly related to the dipole nature of the earth magnetic field for the E and W peaking, and to the GC for the EW 
asymmetry. It is also seen that the EW asymmetry decreases as expected with the increasing energy. 

Eor intermediate geomagnetic latitudes ffigure IT^ . the West peak appears dampen and wider compared to low 
latitudes, while the East peak is even fainter and disappears in the high energy bins. The two-peaks pattern reappears 
with nearly equal heights at high latitude (figure^)). This latter symmetry was expected since the GC disappears in 
the polar region. 

Eor the azimuth angle distributions of the flavor ratios, it can be noted that the z/e/z^e ratio is strongly EW 
asymmetric with a broad enhancement on the West side at low and intermediate latitudes. The opposite trend is 
observed for the ratio, while the (z/^ + + ^e) ratio is found to be almost structureless at all latitudes. 

The reason for these features is similar to that of zenithal angle dependence. The latitude dependence of the 
zenith and azimuth angle distributions reflects the characteristic of the geomagnetic field and reminds us that an 
appropriate treatment of the geomagnetic field is required to account for the upward going muon event data and 
for the east-west asymmetry. Besides, detailed simulation for detectors at different locations is needed to provide a 
sensitive investigation of the geomagnetic dependence of the measurements. 


G. Cosmic Ray parent proton energy distributions 


The incident CR proton energy distribution inducing the neutrino flux are shown on figure for neutrino energies 
around 0.3, 1, 3, and 10 GeV. The distributions are given for + z7g and z^^ + z^. They are averaged over the whole 
detection sphere and Att solid angle, and presented as cumulative probabilities for neutrinos within the energy bins to 
be produced by protons with an energy below the value given in abscissa, and energy distributions normalized to 1. 

At variance with the results of |T3, the mean energy of the primary proton flux producing 0.3 GeV (z^e + ^e) 
obtained in these calculations, appears to be hardly smaller than the for {y^ + in the same energy bin. The 
difference becomes larger with the increasing energy however. The mean values of the proton distributions of the 
figure are given in table ITTTl 

Table ITVl gives for comparison, the mean energy of the {ue + ^e) and neutrino spectrum produced by 

protons within the indicated energy bins. The mean neutrino energy value is found to be smaller for (z^g + ^e) than for 
(^u + for the higher energy bins of the incident protons. This is consistent with the fact that in the decay chain 
TT^ ^ /i^ -hz/g(z/g) + z/^(z/^), the z/g(z/g) particle shares the energy of decaying /i^ which momentum 

is equal to that of the first in tt^ center of mass frame with and another z^^(z^^). The energy-momentum 
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FIG. 10: Zenith angle distribution of the neutrino flux (top) and of the flavor ratios (bottom), in the indicated energy bins 
and for the indicated flavor, for low geomagnetic latitudes {Oiat < 0.5 rad), and averaged over An solid angle. Here, the zenith 
angle is defined such that cos 0 = 1 corresponds to the downward direction. 


E^{GeV) 

0 . 1 - 0.2 

0 . 28 - 0.32 

0 . 8 - 1.2 

2 . 8 - 3.2 

8 . 0-12 

20-30 

0 . 1-30 

< Ep >ue+ue 

30.1 

37.1 

57 

102 

220 

405 

39 

< Ep > vp^-\-i>p 

31 

39 

64 

120 

271 

498 

43 


TABLE III: Mean incident kinetic energy in GeV of the GR protons producing neutrinos in the energy bins Ej^. 
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Zenithoi ongie dtstnbutiorts at medium geomogoetic fotUude 
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FIG. 11: Same as figure IT()1 for the geomagnetic latitudes 


0.5 < Oiat < 1. 


Bp (GeV) 

0.2-10 

10-30 

30-50 

50-200 

200-500 

500-2000 

0.20-2000 

< Bt^e + i'e ^ 

0.25 

0.40 

0.52 

0.67 

0.93 

1.14 

0.42 


0.25 

0.40 

0.56 

0.81 

1.36 

2.19 

0.48 


TABLE IV: Mean neutrino energies <Eu> produced by CR protons within kinetic energy bins Ep. Neutrinos with energy 
below 0.05 GeV are not included in the average. 
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ZenUhoi ongle distnbuttons ot high geonnagneiic lotUude 
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FIG. 12: Same as figure E3 for the geomagnetic latitudes 


1 < Oiat < 7r/2. 


conservation implies that the z^e(^e) energy is on the average smaller than the mean energy of the two 

On figure El it is seen that low energy primary protons contributes more to {ue + ^e) than to the same energy 
+ ^fi) production. This is easy to understand if we note that higher energy muons which are produced by higher 
energy primaries on the average, do not contribute to electron neutrino, but accompany with them, a muon neutrino 
has been created in pion decay. Even at low energy, this process can happen near the detection surface, and so the 
mean primary proton energy for muon neutrino production is raised up. 
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AzJmothol ongle distrlbuttons ct low geomognetic lotUude 
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FIG. 13: Azimuth angle distribution of the neutrino flux (top) and of the flavor ratios (bottom), in the indicated energy bins 
and for the indicated flavor, for low geomagnetic latitudes {Oiat < 0.5 rad), and averaged over An solid angle. The angle 
convention is 0 = 0^, 90^ for geomagnetic south and east, respectively. 


H. Neutrino production altitudes 

The distribution of the neutrino production altitude is plotted in figure El where about 90% neutrinos appear to be 
produced between 5 and 40 km altitude, with the most probable production altitude being within the 15-20 km range. 
These numbers are consistent with the previous 1-dimensional analysis [^. Below a few kilometers, the production 
distribution tends to flatten, governed by the absorption rate, while beyond 40 km, it drops following the atmospheric 
density profile. 
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Table M lists the average production altitudes corresponding to different neutrino energy bins. The average electron 
neutrino production altitude appears to be a little below that of muon neutrinos, because muons can ffy over a certain 
distance - on the km scale for the energies considered here - before decaying into electrons. 

The distribution of the altitude where the first interaction of the incidents CRs with atmosphere takes place was 
found to have a mean value of about 26 km on average, 23 km for nearly vertical (0.9 < cos^ < 1.0 at the first collision 
point) and 38 km for nearly horizontal (0.0 < cos^ < 0.1 at the first collision point) CRs, with no other selection cut 
than a produced neutrino reaching the detection sphere eventually. 

FigurelT^ shows the distribution of the averaged altitude at which the N-th collision occurs in the neutrino producing 
reaction chain as a function of the number N of the collision. 
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FIG. 14: Same as figure IT^ for intermediate geomagnetic latitudes (0.5 < Out < !)• 
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Azimuthol angle distfibuttons ot high geomognetic totUude 
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Neutrino flovor rotio ot high geomognetic lotitude 


o 

o 

X 

i 2 

1 

0,1<E,<0.3GeV 

. , , , 1.1 , . 

0.3<E,<1.0 GeV 

. 1 , , , , 1 , , . . 1 . , 

1.0<E,<3.1 GeV 

, 1 , , . , 1 . , , , 1 , , 

3.1<E,<10GeV 

, 1 , . . , 1 , . . , 1 . . 

g 

g 

51.25 

C 

0.1<E,<0.3 GeV 

; 0.3<E,<1.0GeV 

1.0<E,<3.1 GcV 

—-n — 

3.1<E,<10GeV 

rr 



— . _— -- — 


- r+^ 

0.75 

... , 1 , , 

, 1 ,. . ..1. . 

■..... 1,,,, 1,, 

- ^ 

■..1. , 

*5 3 

o 

"0.1<E^<0.3 GeV 

"0.3<E,<1.0GeV 

"1-0<E^<3.1 GeV 

"il<E,<10GeV 

X 

3 

u: 2.5 





2 

1.5 

~ 1 M 1 1 M , h 1 1 1 1 1 , 

~ ^ 1 i I i 11,1 i i 1111 1 

~ i M i 1 i i I I M I I I i ■ 

" ^ _ 

, 1 111.11111111 i 111 


' I I I i I 1 I I i I h I I I I I i I M I i i i I I i I i I I I I I I I I i I I i 1 i i I 1 M I I 1 i ■ I , I I I 1.1 I I I 1 I I I i I I 1 

0 100 200 300 0 100 200 300 0 100 200 300 0 100 200 300 


^ P 


FIG. 15: Same as figure IT^ for high geomagnetic latitudes (1. < Oiat < 7r/2). 


It can be observed that, on the average, the second collision takes place around 19 km, and the third collision 
around 15 km. In addition, below 4 km altitude, there is almost no meson producing collision because the initial 
CR energy has been damped through the cascade and most particles produced below this altitude are below the pion 
production threshold. 
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Primary pfoton energy (GpV) 


FIG. 16: The cumulative (top) and normalized (bottom) primary proton energy distribution for z^e + Fe (solid lines) and 
(dashed lines) production in 5 energy bins 0.28-0.32, 0.8-1.2, 2.8-3.2, 8.0-12.0 and 20.0-30.0 GeV from the left to the right 
respectively. 


E^iGeV) 

0.1-0.2 

0.28-0.32 

0.8-1.2 

2.8-3.2 

8.0-12 

20-30 

< Altjyg-i-i/g > 

18. 

18.5 

18.3 

17.1 

16.1 

15.3 

< Alti/^+i?^ > 

18.1 

19 

19.6 

19.6 

19.3 

19.1 


TABLE V: Mean neutrino production altitudes (km) corresponding to different neutrino energy bins. 


I. Neutrino flying distance 

Figure shows the distribution of the neutrino flying distance between production and detection point. The 
distributions of the mean values as a function of the cosine of the zenith angle are plotted on figure (201 for four 
neutrino energy bins and for the four neutrino flavors. These distributions are very similar to those obtained in M 
The degeneracy of the curves shown in figure (201 for the 4 flavors can be easily understood by looking at figure ITTI 
and table El the scales of the peak production altitude and the difference of the production altitude for different 
flavors are so small in comparison with the earth diameter that it can be negligible. 


J. 3-Dimensional effects 

The normalized distributions of angles between the directions of the incident CRs at the first collision and of the 
neutrinos in the end of the reaction chain, averaged over the whole detection sphere and 4 tt solid angle, are plotted 
on figure [?T1 for three energy bins. 

The mean values corresponding to the different energy bins are listed in table ED These values indicate that the 
3-dimensional effects are not important for GeV neutrinos, which is consistent with the general expectation and the 
result reported in 

The mean angles between the momenta of: 1)- the incident CR proton (at the first collision) and the pion-producing 
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FIG. 17: Cumulative (left) and normalized (right) and production altitude distributions, for different neutrino 

energy bins, averaged over the whole detection sphere and Att solid angle. 



FIG. 18: Average collision altitude of the iVth collision versus collision number N. 
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FIG. 19: Cumulative (left) and normalized (right) z^e, and Hying distance distributions, for different neutrino energy 

bins, averaged over the whole detection sphere and Htt solid angle. 


secondary proton (at the nth collision in the cascade), 2)- The incident proton (called pion parent) and the produced 
pion in the production collision (production angle), and 3)- The muon parent (tt or K) and the produced muon, are 
given in tables ivm and ivTm for different neutrino energy bins. From these numbers it appears that large deviations 
of the particle directions occur at the pion production level and only at low energy. 


(GeV) 

0.1-0.2 

0.28-0.32 

0.8-1.2 

2.8-3.2 

8.0-12 

20-30 

0.1-30 

<^CR-(ue,ue)> 

28.3 

16.7 

7.3 

3.2 

1.4 

0.69 

20.2 


27.7 

15.7 

6.6 

2.7 

1 

0.46 

19.1 


TABLE VI: Mean neutrino production angle (deg) relative to the CR particle direction at the first collision, for different 
neutrino energy bins. 


(GeV) 

0.1-0.2 

0.28-0.32 

0.8-1.2 

2.8-3.2 

8.0-12 

20-30 

0.1-30 

^^CH-7r±par.> 

4.98 

2.90 

1.22 

0.57 

0.32 

0.26 

3.56 

^?7r±par.-7r±^ 

24.28 

14.09 

5.88 

2.33 

0.85 

0.40 

17.28 


3.73 

1.95 

0.72 

0.30 

0.26 

0.25 

2.55 


TABLE VIL Mean angles between muons, pions, and their parent particle (p,n,K), and between the parent particle and the 
incident CR, for different energy bins of z^e(i^e) (deg). 


The deflection of the charged particles by the geomagnetic field has also been examined. The largest deflection is 
undergone by the primary proton between the generation point and the first collision point (which depends on the 
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FIG. 20: Mean flying distance versus the zenithal angle distribution of z^e, and in the indicated energy bins (GeV), 

averaged over the whole detection sphere. 


+ (GeV) 

0.1-0.2 

0.28-0.32 

0.8-1.2 

2.8-3.2 

8.0-12 

20-30 

0.1-30 

^^CH-7r±par.> 

5.33 

2.91 

1.18 

0.54 

0.30 

0.26 

3.89 

^^7r±par.-7r±^ 

23.96 

13.41 

5.36 

1.99 

0.66 

0.32 

16.48 


4.00 

2.12 

0.79 

0.30 

0.26 

0.25 

2.73 


TABLE VIII: Same as table lyTTI for different energy bins of . 


generation radius). The pion deflection angle is negligible in account of its short lifetime, while for muons it is less 
than 4°. 

Hence, it is confirmed that the 1-dimensional approximation used for GeV neutrinos in the previous calcula¬ 
tions / simulations was an acceptable approximation. 


V. NEUTRINO FLUX AROUND THE SUPER-KAMIOKANDE DETECTOR 


In this section, the features of the atmospheric neutrino flux at the geographical position of the Super-Kamiokande 
(SuperK) detector are investigated. The detection sphere was defined at 0.372 km altitude, and within (36'^25'33" ± 
l.b^)N and (137^18'37" ± lb^)E latitude and longitude bins respectively. 


A. Flux and flavor ratio 

The calculated energy distributions of the atmospheric neutrino flux and the v/v flux ratios around the SuperK 
detector, averaged over Air solid angle, are compared with the results of Honda et al. 1^ ^ on figure [221 for the 
various flavors. The flux obtained in the present work appear to be significantly smaller than in the 1-Dimensional 
calculations for low neutrino energies. 
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FIG. 21: Normalized distributions of neutrino production angle relative to the direction of the incident CR proton at the first 
collision, for three neutrino energy bins (GeV) : 0.28-0.32, 0.8-1.2, 2.8-3.2. 


The values of the simulated 0.1-20 GeV neutrino flux for the location of the Super-Kamiokande experiment are 
tabulated in table ITYl These spectra can be fit by means of the function form: 


fiE,) = Cl (1.0 + C 2 Exp( -C3 E, )) (1) 

very nicely, where f{Ey) is the flux, Ey the neutrino energy, and ci,C 2 ,C 3 and C 4 are fitting parameters given in 
table 13 

The functional form is inspired from carrying some features of the hadroproduction of pious and kaons. 

For the flavor ratios, our (z/^ + ^^)/(^e + ^e) is similar to that reported in gm, but z^e/^e differ evidently from 
all the previous 1 -dimensional calculation summarized in figure 12 of 0 . 

From figures Q] and El it can be observed that the positive muon flux is overestimated in our simulation, which 
subsequently raises the z^e/^e ratio. This observation can explain partially the large difference observed. 

However, as explained at the end of section imi it is believed that the tendency of decreasing z/g/z^e with the 
increasing energy should be expected because of the constraint on the j' k~ [K^ jK~) ratio imposed by the relevant 
data. 


B. Zenith angle distribution 

The zenith angle distribution of the flux and of the flavor ratios are shown on figure [ 2 SI where the enhancement of 
the flux in horizontal directions discussed in section rrvu is observed, in qualitative agreement with M- In the low 
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K (GeV) 


FLUX 

(m-^s-isr-^GeV^ 






0.100-0.124 

2948 ± 84 

1780 ± 59 

5229 ± 103 

4915 ± 99 

0.124-0.153 

2078 ± 53 

1324 ± 45 

3677 ± 68 

3629 ± 80 

0.153-0.189 

1595 ± 50 

892 ± 29 

2771 ± 63 

2569 ± 48 

0.189-0.233 

1023 ± 30 

702 ± 29 

1875 ± 40 

1853 ± 40 

0.233-0.289 

774 ± 28 

485 ± 19 

1260 ± 27 

1266 ± 27 

0.289-0.357 

503 ± 15 

335 ± 10 

909 ± 20 

901 ± 21 

0.357-0.441 

325 ± 12 

238 ± 8 

649 ± 17 

610 ± 14 

0.441-0.545 

218 ± 8 

162 ± 7 

450 ± 15 

389 ± 9 

0.545-0.674 

140 ± 5 

113 ± 5 

291 ± 9 

268 ± 7 

0.674-0.833 

99 ± 6 

70 ± 3 

185 ± 5 

183 ± 5 

0.833-1.029 

54 ± 2 

42 ± 2 

125 ± 5 

114 ± 4 

1.029-1.272 

33 ± 1 

26 ± 1 

72 ± 2 

70 ± 2 

1.272-1.572 

27 ± 2 

18 ± 1 

46 ± 2 

46 ± 2 

1.572-1.943 

13.1 ± 0.8 

13 ± 2 

29 ± 1 

26 ± 1 

1.943-2.402 

8.8 ± 0.6 

5.7 ± 0.5 

20 ± 2 

17.4 ± 0.8 

2.402-2.969 

4.1 ± 0.3 

3.9 ± 0.4 

9.7 ± 0.4 

10.1 ± 0.4 

2.969-3.670 

2.4 ± 0.2 

1.8 ± 0.2 

6.4 ± 0.4 

5.5 ± 0.3 

3.670-4.537 

1.2 ± 0.2 

1.1 ± 0.1 

3.8 ± 0.2 

3.1 ± 0.2 

4.537-5.608 

0.79 ± 0.08 

0.54 ± 0.05 

2.2 ± 0.1 

3.1 ± 0.7 

5.608-6.931 

0.35 ± 0.04 

0.27 ± 0.03 

1.14 ± 0.07 

1.2 ± 0.1 

6.931-8.568 

0.19 ± 0.02 

0.19 ± 0.04 

0.55 ± 0.04 

0.72 ± 0.06 

8.568-10.59 

0.12 ± 0.03 

0.22 ± 0.07 

0.47 ± 0.03 

0.32 ± 0.02 

10.59-13.09 

0.033 ± 0.005 

0.036 ± 0.005 

0.17 ± 0.01 

0.158 ± 0.009 

13.09-16.18 

0.035 ± 0.006 

0.022 ± 0.004 

0.103 ± 0.007 

0.086 ± 0.007 

16.18-20.00 

0.013 ± 0.003 

0.009 ± 0.002 

0.046 ± 0.004 

0.040 ± 0.003 

TABLE IX: Neutrino flux calculated for Super-Kamiokande experiment location. 

Particle 

Cl 

C2 

C3 

C4 

Ve 

89.92 

-1.043 

0.776 

-2.99 

Ve 

60.47 

-1.056 

0.935 

-2.88 


198.4 

-1.029 

0.744 

-2.82 


201.3 

-1.031 

0.674 

-2.88 


TABLE X: Parameters used to fit the simulated neutrino spectra around the Super-Kamiokande detector. 


energy bin, the maximum of this enhancement is about twice the downward or upward going flux out of the horizontal 
plane in the low energy bin. Again, it can be seen on the figure (right panel), that the flavor ratio + ^e) 

also depends on the zenith angle, especially for high energy neutrinos, for the same reason as given previously. 


C. Azimuthal angle distribution 

The azimuth angle distribution of the flux and of the flavor ratios are shown on figure for the four neutrino 
species. The large angular binning used in the figure was dictated by the low counting statistics obtained, even for 
very long computer time used for running the simulation program. The East-West effect discussed previously in 
section HV FI is clearly observed here, more prominent for low energy neutrinos. The EW asymmetry discussed below 
has been estimated on the basis of these distributions. 


D. Estimate of the East-West asymmetry 

In addition to the atmospheric muon flux, the EW asymmetry measurements provide an important global test of 
the reliability of the overall approach. Since muons are charged particles, having finite lifetime, and losing energy 
during their flying in the atmosphere, while neutrinos have none of these features, these differences make the relation 
between muon and neutrino flux a complicated issue 0- However, the EW asymmetry is independent on the physics 
beyond the standard model and is not expected to depend on the neutrino mass or oscillations. A reliable simulation 
should thus be able to reproduce the experimental EW asymmetry data quantitatively. 
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FIG. 22: Simulated atmospheric neutrino spectra and flavor ratio around Super-Kamiokande detector. 
For the asymmetry parameter Aew defined as: 


_ Ne — Nw 
Ne + Nw 


Ne{Nw) being the number of eastward (westward) lepton events, the SuperK Collaboration reported the following 
measured values: 


= 0.21 ± 0.04 = 0.08 ± 0.04 

for e-like and /i-like events respectively, for a selection of single-ring events with momentum between 0.4 — 3 GeV and 
—0.5 < cos Tzenith < 0-5, with Tzenith being the zenith angle. 

In order to obtain the Aew value precisely, the 3-Dimension differential neutrino reaction cross section 
needs to be known. The asymmetry Aew however could be estimated by using the following procedures: 

1. Evaluate the energy of the neutrinos which can induce lepton events over the same range as selected by superK, 
i.e., 0.4 to 3 GeV. The leptons with momentum pi ^ 100 MeWc, appear to carry 65% of the incident neutrino 
energy, this proportion increases to ^ 85% for pi = I GeV/c It can be assumed that the fraction keeps on 
increasing and asymptotically tends to 100%. Based on this supposition, the corresponding neutrino energies of 
about 0.55 and 3.1 GeV associated to the 0.4 and 3 GeV lepton events are obtained respectively. 

2. Estimate the production angle of the neutrino induced lepton. This can be done following [^, by noting that 
the interested neutrino energy range is above 0.5 GeV within which there is no large difference between the 
production angles of Ue and events. Another way is by using the values 55^ at pi = 0.4 GeV/c and 20^ at 
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FIG. 23: Zenith angle distribution of the calculated neutrino flux, and of the corresponding flavor ratios around the Super- 
Kamiokande detector. 

Pi = 1.5 GeV/c as reported in and fit them as a function of the neutrino energy, assuming that the angle 
smoothly tends to zero as expected from basic kinematics. 

3. For each neutrino with energy between 0.55 and 3.1 GeV, lepton events are randomly generated with the 
scattering angle relative to the neutrino direction estimated as above, and with a uniform distribution for 
the azimuth angle. After weighting by the energy dependent cross section mis EH and application of a 
normalization procedure, the events with | cos{0i) \< 0.5 are selected to calculate the asymmetry parameter. 

Using the above procedure, the following values are obtained for the EW asymmetry parameters for e-like and 
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Neutrino flovor rotio oround Super-Komiokonde 



FIG. 24: Simulated azimuth augle distributiou of the ueutriuo flux aud of the flavor ratios arouud the SuperK detector. 


/i-like events around the Super-Kamiokande detector: 

= 0.12 ± 0.03 = 0.13 ± 0.02 

In contrast, the EW asymmetries obtained for the two neutrino flux with the same cut over the zenithal angle and 
within the same energy bin of [0.55,3.1] GeV, are: 

= 0.17 ± 0.03 = 0-22 ± 0.02 

which shows that the predicted neutrino asymmetry is largely washed out by the angular distribution of the neutrino 
induced lepton production process. 
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The larger EW flux asymmetry obtained for than for z^e(^e) can be qualitatively understood: The production 

of z^e(^e) takes place in the second step of the pion decay chain. The decay kinematics should thus wash out the EW 
asymmetry of neutrinos more than for muons which are produced dominantly in the first stage of the decay sequence. 

In comparison with the data and other simulation result HH , the present calculations departs further from the 
data. It can be noted that, the present overestimated ratio of tt^/tt” may affect this prediction. More works is still 
needed to clear this problem. 


VI. NEUTRINO FLUX AROUND THE SOUDAN EXPERIMENT 

The same calculations for the simulated atmospheric neutrino flux around the SOUDAN detector as for superK 
experiment in the previous section, are reported in this section. The normalization area used in the calculations 
corresponds to the geographical latitude and longitude bin of (48^ ± 5^)N and (98^ ± 10^)IT respectively. The 
altitude is taken the same as for Super-K, as we have seen from figure El the error resulting from the altitude 
difference between Super-Kamiokande and Soudan detector should be negligible in evaluating the flux and relevant 
distributions. 


A. Flux and flavor ratio 


The energy spectra and energy dependence of flavor ratios, averaged over An solid angle are shown in figure E) Due 
to the relatively high geomagnetic latitude of the the Soudan detector location, and hence to the lower geomagnetic 
cutoff, the calculated neutrino flux are much higher than found for Super-K. 


(GeV) 


FLUX 

(m-^s-^sr-iGeV^) 


ye 

ye 

y^ 

y^ 

0.100-0.124 

6222 ± 240 

3344 ± 163 

10967 ± 313 

9453 ± 257 

0.124-0.153 

4338 ± 175 

2283 ± 97 

7154 ± 194 

6604 ± 188 

0.153-0.189 

2698 ± 100 

1674 ± 83 

5078 ± 159 

4872 ± 154 

0.189-0.233 

1738 ± 71 

997 ± 49 

3344 ± 110 

3254 ± no 

0.233-0.289 

1200 ± 47 

860 ± 62 

2062 ± 60 

2124 ± 78 

0.289-0.357 

771 ± 32 

623 ± 45 

1396 ± 53 

1397 ± 52 

0.357-0.441 

488 ± 26 

361 ± 19 

936 ± 32 

876 ± 36 

0.441-0.545 

333 ± 20 

214 ± 13 

628 ± 28 

595 ± 31 

0.545-0.674 

212 ± 15 

163 ± 21 

359 ± 15 

352 ± 15 

0.674-0.833 

122 ± 9 

92 ± 7 

231 ± 10 

234 ± 12 

0.833-1.029 

73 ± 6 

60 ± 4 

145 ± 6 

138 ± 7 

1.029-1.272 

48 ± 3 

34 ± 3 

80 ± 4 

92 ± 6 

1.272-1.572 

28 ± 3 

19 ± 2 

52 ± 3 

49 ± 3 

1.572-1.943 

13 ± 1 

10 ± 2 

30 ± 2 

32 ± 4 

1.943-2.402 

9 ± 1 

6.1 ± 0.7 

19 ± 1 

17 ± 1 

2.402-2.969 

4.4 ± 0.5 

3.5 ± 0.5 

9.9 ± 0.7 

11.1 ± 0.9 

2.969-3.670 

3.4 ± 0.5 

1.5 ± 0.2 

5.5 ± 0.5 

6.4 ± 0.6 

3.670-4.537 

1.6 ± 0.3 

0.9 ± 0.1 

3.6 ± 0.3 

3.3 ± 0.4 

4.537-5.608 

0.53 ± 0.08 

0.49 ± 0.08 

2.0 ± 0.2 

1.5 ± 0.2 

5.608-6.931 

0.44 ± 0.07 

0.43 ± 0.09 

1.1 ± 0.1 

0.9 ± 0.1 

6.931-8.568 

0.16 ± 0.03 

0.12 ± 0.02 

0.58 ± 0.05 

0.6 ± 0.1 

8.568-10.59 

0.10 ± 0.02 

0.08 ± 0.02 

0.32 ± 0.03 

0.29 ± 0.03 

10.59-13.09 

0.031 ± 0.007 

0.04 ± 0.01 

0.20 ± 0.02 

0.16 ± 0.02 

13.09-16.18 

0.017 ± 0.004 

0.012 ± 0.004 

0.10 ± 0.01 

0.09 ± 0.01 

16.18-20.00 

0.016 ± 0.005 

0.022 ± 0.009 

0.053 ± 0.006 

0.037 ± 0.005 


TABLE XI: Neutrino flux calculated for Soudan experiment location. 


Also, the 0.1-20.0 GeV neutrino energy spectra simulated for the Soudan experiment site are tabulated in table IXll 
They can be fit with the same functional form (relation^. The coefficients are listed in table nrm 
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Particle 

Cl 

C2 

C3 

C4 



105.0 

- 1.034 

0.869 

- 3.12 


63.9 

- 1.069 

1.235 

- 2.99 


143.0 

- 0.970 

1.859 

- 2.67 


362.2 

- 1.021 

0.416 

- 3.16 


TABLE XII: Fitting parameters for simulated neutrino spectra around Soudan detector. 





FIG. 25: Simulated atmospheric neutrino spectra and flavor ratios around the SOUDAN detector. 


B. Zenith angle distributions 

The calculated zenithal angle distributions are shown in figure [2^ for the same energy bins as previously, where 
the high geomagnetic latitude feature of such distribution illustrated in Section HV El is evident. 


C. Azimuth angle distributions 

Figure E3 shows the azimuth angle distributions of the atmospheric neutrino flux around the SOUDAN detector, 
based on which the EW asymmetry parameter for neutrino flux and lepton events in the detector are also estimated 
along the same lines as for the superK experiment. 

The EW asymmetry for the lepton events is predicted as 

A^e-w" = 0-006 ± 0.037 A^'e-w" = 0-026 ± 0.024 
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Zenithol ongle distribution oround Soudan detector 



Neutrino flovor ratio oround Soudon detector 



FIG. 26: Zenith angle distributions of the simulated neutrino flux and flavor ratios around the SOUDAN detector. 


while for the neutrino flux the following values are obtained: 

= 0.068 ± 0.043 = 0.071 ± 0.025 

with the same selection criteria as used for SuperK in the previous section. It would therefore be of great interest to 
have some experimental values for this experiment, to provide a further test the present approach, and a step forward 
to a better global understanding of the problem. 
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Azimuthol angle distributton oround Soudon detector 



Neutnno flovor ratio oround Soudon detector 



FIG. 27: Azimuth angle distributions of the simulated neutrino flux and of the flavor ratios around the SOUDAN detector. 


VII. CONCLUSION AND DISCUSSION 

In summary, the atmospheric neutrino flux has been simulated by means of an event generator dedicated to the 
general Cosmic Ray - Atmosphere interactions, developed to account for the AMSOl results, originally. The successful 
calculations reported previously on the hadron and lepton flux close to earth, strongly support the method and the 
models used. 

In this work, the CAPRICE and HEAT muon flux data measured at various altitudes in the atmosphere, from about 
sea level to 38km have been successfully reproduced, providing further grounds to the correctness of the approach. 

A significantly lower absolute value of the flux than found in 1-dimensional calculations 11113 in the low energy 
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range has been reported here. Also, the Vejve ratio is largely different from all the previous calculations. In addition, 
the present results are in agreement with similar results reported in 

Detailed features of the neutrino ffux: zenithal and azimuthal angle distribution and ffavor ratios at different 
latitudes, corresponding energy distributions, production altitude and production angle distributions, have been 
calculated and discussed. The components originating from Cosmic Rays and from the atmospheric cascade have 
been distinguished and discussed separately. All the observed features of the calculated ffux could be traced back to 
the primary spectra, 7r+/7r“(i^+/i^“) ratio, muon kinetics, geomagnetic effect, and geometry. 

Specific 3-dimensional effect seen in the simulation results are very small, which indicates that the 1-dimensional 
approximation provides reliable results for GeV neutrino ffux calculations. 

The calculations fail to reproduce quantitatively the East-West asymmetry of the atmospheric neutrino-induced 
events in the Super-K detector. More investigations are needed to understand this disagreement. 

A complete comparison of the parametrizations used to describe the various hadronic production channels with the 
available data will be reported in details later. 
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